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ABSTRACT 

We argue that quiescent WZ Sge is a rapidly spinning magnetic rotator in which most 
of the matter transferee! from the secondary is ejected from the system. Assuming that 
the observed 27.87 s oscillation period is due to the spinning white dwarf we propose 
that the other observed principal period of 28.96 s is a beat due to reprocessing of the 
rotating white dwarf beam on plasma blobs in Keplerian rotation at the outer disc 
rim. The weaker, transient, 29.69 s period is identified as a beat with the Keplerian 
period of the magnetosphere. WZ Sge evolves through a cycle of spin-up and spin- 
down phases. During the spin-down phase it is a DQ Her star, during the spin-up 
phase it should be a ER UMa star. 

Key words: stars: individual: WZ Sge - novae, cataclysmic variables - stars: oscil- 
lations - stars: rotation - stars: magnetic fields 



1 INTRODUCTION 

WZ Sagittae is a remarkable dwarf nova binary system. 
Dwarf novae are a subclass of cataclysmic variables (CVs) 
which show more or less regular outbursts. In CVs, a late- 
type, Roche-lobe filling secondary star loses mass which, in 
general, is accreted by a white dwarf primary. Usually, when 
the white dwarf is not too strongly magnetized the accreting 
matter forms an accretion disc. It is well established that 
such discs are the sites of dwarf nova outbursts. WZ Sge 
outbursts are very rare (the recurrence time is around 30 
years) and they are only of the 'superoutburst' type, i.e. 
they are long (~ 30 days), high amplitude (~ 7 mag) out- 
bursts during which a 'superhump' in the optical light-curve 
is observed at a period slightly longer than the orbital one. 
WZ Sge is the prototype of a class of systems with similar 
outburst properties (sometimes also called 'TOADs'). The 
orbital period of WZ Sge is 81.63 minutes, close to the min- 
imum period (~ 80 min) below which no (non-degenerate) 
CV has been observed. 

It is generally believed that dwarf nova outbursts re- 
sult from a thermal-viscous instability present in accretion 
discs at temperatures corresponding to hydrogen partial ion- 
ization. In the standard version of the disc instability model 
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(DIM; see Hameury et al. 1998 for the most recent version of 
this model, Cannizzo 1993 and Lasota & Hameury 1998 for 
reviews) one assumes that mass-transfer from the secondary 
is constant prior to, and during the outburst. Properties of 
outbursts depend on the viscosity mechanism which trans- 
ports angular momentum and heats the disc. In the DIM 
the kinematic viscosity coefficient is taken as v = ac s H, 
where c s is the (adiabatic) speed of sound, H the disc semi- 
thickness, and the 'viscosity coefficient' a < 1 (Shakura & 
Sunyaev 1973). The DIM requires different values of a's in 
outburst and in quiescence (Smak 1984b). The quiescent 
value for the great majority of systems should be ~ 0.01 
(Livio & Spruit 1991). Smak (1993) showed, however, that 
the outburst cycle of WZ Sge cannot be described by such a 
version of the DIM (see also Osaki 1996). The main reason 
is that for a ~ 0.01 the disc is able to accumulate during 
quiescence only ~ 10 21 g whereas during outbursts ~ 10 24 
g is accreted by the white dwarf. Assuming that a ~ 10 -5 
solves this inconsistency and the resulting very low viscos- 
ity may also account for the long recurrence time which is 
difficult to reproduce with the 'standard' a values. The dif- 
ficulty of this version of the model is that one must find an 
explanation for the very low value of the viscosity coefficient 
in this particular system. The low value of the mass-transfer 
rate (~ 10 15 g s _1 , invoked as an explanation; Osaki 1996) 
is not much lower than in other systems with similar orbital 
periods so it is doubtful that it can explain a difference in 
a by three orders of magnitude. In any case the relation be- 
tween accretion rate and viscosity, if any, is unknown (see 
Gammie & Menou 1998). 

A different modification of the standard DIM was sug- 
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gested by Lasota, Hameury & Hure (1995) and Hameury, 
Lasota & Hure (1997). They assume that the value of the 
a-parameter in WZ Sge is not different from the 'standard' 
one. They noticed that at low mass transfer rates outer re- 
gions of accretion discs in CVs are cold enough to be stable 
with respect to the thermal instability believed to be re- 
sponsible for dwarf-nova outbursts. They proposed that the 
WZ Sge accretion disc does not extend down to the white 
dwarf's surface but is truncated at a radius corresponding 
to a stable outer disc. An enhancement of mass transfer 
would bring such a disc into an unstable state and thus trig- 
ger an outburst. One should stress that such an outburst 
would still be due to the thermal- viscous instability, the en- 
hanced mass transfer serves only as a trigger. The recurrence 
time is then related to the characteristic time-scale of mass- 
transfer fluctuations. One can also obtain long recurrence 
times if the truncated disc is marginally unstable (Warner, 
Livio & Tout 1996) but in practice such a disc is indistin- 
guishable from a marginally stable one proposed by Lasota 
et al. (1995). There is still, however, a problem to be solved: 
since a ~ 0.01 there would be only ~ 10 21 g in the disc 
available. Therefore, in this case mass must be added to the 
disc during outburst. There is observational evidence (Smak 
1997; see also Hessman et al. 1984) that irradiation of the 
secondary during outburst increases mass transfer from the 
secondary (see Hameury et al. 1997 for a simple model). 

In Lasota et al. (1995) and Hameury et al. (1997) the 
truncation was assumed to be due either to evaporation 
(Meyer & Meyer-Hofmeister 1994) or to the presence of a 
magnetic field strong enough to disrupt the quiescent disc 
(iS; 10 4 G would be sufficient for a low mass white dwarf). 
The second hypothesis seemed to be favoured by the pres- 
ence, before the last, 1978 outburst, of a 27.87 s coherent 
oscillation which could be attributed to the rotation of an ac- 
creting white dwarf (Patterson 1980, hereafter P80). There 
were two problems with this interpretation. First, in ad- 
dition to the 27.87 s period several other (some of them 
transient) 'satellite period' oscillations have been observed 
(Robinson, Nather & Patterson 1978; P80). The interpre- 
tation of these periods in terms of beat frequencies result- 
ing from reprocessing of the white dwarf's pulsed light on 
various features in the accretion flow encountered several 
difficulties (P80) . Second, all the oscillations disappeared 
during the outburst. After the outburst the principal, 27.87 
s, oscillation has been absent for 16 years (although the 
28.96 s was seen when WZ Sge was at twice its pre-outburst 
brightness) . Although it was easy to understand that during 
the outburst the increased accretion rate could suppress the 
magnetosphere, the lack of the principal pulse afterwards 
was difficult to understand and cast a shadow of doubt on 
the presence of a rapidly rotating white dwarf in WZ Sge. 

In 1995 the 27.87 s oscillation reappeared in the com- 
pany of the previously present 28.96 s pulsation (Patterson 
et al. 1998, hereafter P98). Weak, transient 'satellite period' 
pulses are also seen at 28.2 and 29.69 s. The most impor- 
tant news, however, was the detection of a 27.86 ± 0.01 s 
period in the 2-6 keV ASCA energy band (P98) which con- 
firmed the presence of a rotating, magnetized white dwarf 
in WZ Sge. The presence of a rapidly rotating white dwarf 
was confirmed by UV spectral observations (Cheng et al. 
1997). WZSge is therefore a DQ Her star, a CV containing 
a rapidly rotating magnetized white dwarf. 



In a recent article Meyer-Hofmeister, Meyer & Liu 
(1998) (see also Mineshige et al. 1998) argue that the in- 
ner disc 'hole' in WZ Sge is due to evaporation and not to 
the presence of a magnetic field. These authors, however, 
underestimate by a significant factor the coherence of the 
observed pulsations and their model seems to be contra- 
dicted by observations. For example, their model predicts 
the presence of superhumps well before the beginning of the 
superoutburst whereas in WZ S ge this feature appeare d only 
10 - 12 day after the maximum ( Patterson et al., 1981 ). Also 
the remarkable constancy of the X-ray luminosity after the 
1978 outburst is in contradiction with their model. 

In this article we argue that Patterson's (P80, also P98) 
hypothesis that WZ Sge contains a rapid, oblique magnetic 
rotator is consistent with most of the observed properties 
of this system. We show that sideband oscillations observed 
in WZ Sge can be interpreted as resulting from reprocessing 
of the white dwarf's pulsed light on various features of the 
transfered flow. In Section 2 we present the observed proper- 
ties of WZ Sge and identify various oscillations with orbital 
sidebands. In Section 3 we discuss possible configurations 
of WZ Sge as a system containing a magnetic rotator and 
show that only one is consistent with the system properties 
at outburst. In Section 4 we discuss the accretion history of 
WZ Sge and its relation to other CVs. We summarize our 
model in Section 5. 



2 WZ Sge AS A MAGNETIC ROTATOR 

If WZ Sge is an oblique magnetic rotator one should, in prin- 
ciple, observe in its optical light-curve the fundamental spin 
frequency and some sidebands resulting from the beat be- 
tween this and the orbital frequencies (Warner 1986). In 
what follows we will assume that all characteristic periods 
observed in WZ Sge have a rotational origin. The possibil- 
ity that these oscillations (or some of them, Warner 1995a) 
are due to white dwarf pulsations cannot be excluded (see 
Wood 1999 for recent arguments in favour of this model). 
Such a model is, however, hard to test because it does not 
make any predictions (P98). In what follows we will try to 
show that all observed periodicities observed in WZ Sge can 
be explained in the framework of a magnetic rotator model. 



2.1 General properties and parameters 

Following Warner (1995a) we classify WZSge as a DQ Her 
star. DQ Her stars are rapidly rotating intermediate polars 
(IPs), i.e. magnetic CVs in which the white dwarf rotation is 
not synchronous with the orbital motion (synchronous sys- 
tems are called 'polars' or 'AM Her stars'). Warner (1995a) 
adds also that DQ Her stars are characterized by the ab- 
sence of hard X-rays in the sense that, contrary to 'usual' 
IPs their X-ray temperature is much lower than 10 keV. In 
WZSge the X-ray temperature is ~ 4.5 keV (P98), so it is 
a hard X-ray emitter, even if its temperature is lower than 
that in 'usual' IPs (see below). 

WZ Sge has been one of the best observed CVs but 
the values of its fundamental parameters are still uncer- 
tain. Smak's (1993) photometric solution gives Mi = 0.45 
(where Mi is the white dwarf mass in solar units) and 
q = MijM\ = 0.13, where Ma is the secondary's mass, 
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Figure 1. Various radius vs. mass relations appropriate for WZ 
Sgc. Shown are the mass-radius (M-R) relation for helium white 
dwarfs (Nauenberg 1972; dashed line), the corotation radius (Rci; 
Eq. 1), the critical radius (-R cr it) above which a disc is marginally 
stable (Eq. 6) for 3 values of a (see text; dash-dotted lines), and 
the Kepler radii at which the beat between the white dwarf spin 
frequency and the Kepler frequency gives rise to the observed 
oscillations at 28.96 s and 29.69 s. The inferred solution from Smak 
(1993) is also indicated. With dotted lines we indicate the white 
dwarf masses as derived by Smak (1993) and Spruit & Rutten 
(1998), i.e. 0.45 M Q and 1.2 M , respectively. 



whereas Spruit & Rutten (1998), who model the 'hot spot' 
at which the mass transfer stream interacts with the outer 
disc regions, get Mi = 1.2 and q = 0.075. As we shall see 
Smak's value is too small if the white dwarf in WZ Sge is 
spinning at 27.87 s, but the Spruit & Rutten (1998) value 
might be too high. 

In Figure 1 we show various mass-radius relations rel- 
evant for WZ Sge. The white dwarf mass-radius relation is 
that of Nauenberg (1972) which is suitable for helium white 
dwarfs (using other mass-radius relations gives very similar 
results, e.g. P98). Also plotted is the corotation radius 

/ GMqMiP^\ ^ g 2/3 , f l/3 /,-, 

Rn = I 4^2 ) =1-5 X 10 PJ Mj/ cm, (1) 

where P 3 is the white dwarf spin period. The corotation 
radius corresponds to a distance at which a free particle in 
circular Keplerian orbit corotates with the white dwarf. In 
particular the white dwarf radius must satisfy Ri < Rn. 
If one assumes that P s = 27.87 s, Figure 1 shows that for 
Smak's Mi = 0.45 the white dwarf would be rotating just 
below the break-up speed so Mi must be larger than 0.45 
(see also P98). 



2.2 Orbital sidebands 

In addition to the 27.87 s oscillation several other optical 
oscillations between 27.8 and 30 s have been observed in 
WZSge (P98, and references therein). In general either the 
27.87 or 28.96 s oscillations are present, sometimes with ad- 



ditional periods up to ~30s (i.e. 28.190, 28.52, 29.69 s). A 
28.2 s oscillation is also seen in the UV (Welsh et al. 1997) 
(strictly speaking the UV oscillation is 28.09 ± 0.12, but we 
will follow P98 in identifying it with the optical 28.2 s (28.19 
s) period). 

We will denote the white dwarf angular rotation fre- 
quency with uj = 2n/P a and the orbital frequency as Q. — 
27r/P or b. As shown by Warner (1986), in intermediate polars 
the rotating white dwarf's beam and its reprocessing on the 
surface layers of the disc, hot spot and/or the secondary pro- 
duces peaks in the power spectra of optical light curves not 
only at uj and uj — fi, but there also exist sidebands at uj+Q 
and UJ—2Q. (Wynn & King 1992 analyzed power spectra of 
X-ray light curves in IPs in the case of discless accretion.) 
For a white dwarf rotation period of 27.87 s and an orbital 
period of 81.63 min (see P98) this gives sidebands at 28.03, 
27.71 and 28.19 s. Patterson (1980) already noted that the 
28.2 s period is consistent with u—2Q. He also noted that 
significant peaks in the power spectra exist in the sidelobes 
close to the main peak at 27.87 s. These periods were, how- 
ever, not discussed by him; visual inspection of his Figure 3b 
shows that the peaks of these sidelobes are at ~27.96 and 
~27.79 s. These are very close to uj— Q and ui+Q, respec- 
tively. We therefore propose that the three optical sideband 
frequencies are present, making stronger the suggestion that 
an appreciable magnetic field exists in WZ Sge, truncating 
the inner parts of an accretion disc and making it similar to 
an IP. It also confirms that the 27.87 s period is the white 
dwarf spin period; checks made by assuming the white dwarf 
spin period is at the other observed frequencies show that 
although the expected sideband frequencies match the ob- 
served frequencies there is the no explanation for the 27.87 s 
period. For example, assuming 28.19 s as the white dwarf 
spin period gives sidebands at 28.03, 28.35, 28.52 s. These 
are all observed in the power spectra of P80, P98. In addi- 
tion uj — 3f2 gives 28.96 s. P s — 28.09 s gives sidebands at 
28.25, 27.93, 28.42 s, which in principle could match the ob- 
served peaks in the power spectrum, as well as uj — 3Q, which 
gives 28.58 s. In both cases, however, the 27.87 s oscillation 
remains unexplained. 

We identify the 28.2 s optical and UV oscillations with 
uj—20. which would result from reprocessing on the hot spot. 
Welsh et al. (1997) identify this period with the white dwarf 
spin period but the presence of the 27.87 s period in X-rays 
does not support this identification. If the 28.2 s oscillation 
results from reprocessing on the hot spot one must then 
explain why the fundamental spin period (27.87 s) is not seen 
in the UV. In order to understand this absence one should, 
presumably, understand the structure of the accretion flow. 
As we shall see below the spin period of WZ Sge implies 
that this system is in an 'ejector' phase, ejecting most of the 
transferred matter. Such a case would be totally different 
from the usual IPs (see e.g. Hellier 1996). No model of such 
a flow seems to exist at present. 

It is also interesting to note that a reprocessing pro- 
jected area might also vary with the 2£l component (Warner 
1995a). This will introduce components at frequencies uj+Q 
and uj— 3fl in the power spectrum. The latter produces a 

3 The 28.19 s oscillation might be composed of signals at 28.14 s 
and 28.24s (P80). 
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peak near 28.35 s, which is visible (at ~ 28.34 s) as a shoul- 
der of the peak near 28.19 s in the 'grand average' power 
spectrum of the observations prior to the outburst end of 
1978 (P80; see also P98). We note that the observed oscil- 
lation at 28.52 s is consistent with oj— 4£1 (P80) but have no 
interpretation for this fact. 

Having explained to our satisfaction various sidebands 
we are still left with the 28.96 s and the 29.69 s periods which 
are not commensurable with the orbital period. The 28.96 s 
signal appears as one of the two main peaks in the power 
spectrum and is generally present. The 29.69 s is weak and 
transitory. 



3 WZ Sge AS A MAGNETIC EJECTOR 
3.1 A discless system ? 

Patterson (1980) suggested that the 28.96 s period could re- 
sult from the beat between the white dwarf spin frequency 
and the Keplerian frequency at the magnetosphere. The re- 
quired Keplerian period is Pk = 733.47 s and could cor- 
respond to the rotation period of plasma blobs at the disc 
inner edge. The problem with such an interpretation is that, 
if true, not much of a disc would be left. The corresponding 
magnetospheric radius would be 



0.45 M., 



J?a/ (28.96) 



1.22 x 10 10 M 1 1/3 cm. 



(2) 



According to Smak (1993), for Mi = 0.45, the outer disc 
radius is R D = (1.07 ± 0.19) x 10 10 cm, whereas Spruit & 
Rutten (1998) obtain, for Mi = 1.2, _R D = (1.75 ± 0.14) x 
10 10 cm so that in both cases we would rather have a ring 
then a disc, or no disc at all, depending on the details of 
plasma interaction with the magnetic field (see e.g. King 
1993). In addition, if the Keplerian radius Rk corresponding 
to the 28.96 s period were the magnetospheric radius, the so 
called 'fastness parameter' (see e.g. Frank, King & Raine 
1992) would be 

Pk 



= 26.3. 



(3) 



Since u s S> 1, matter transferred from the secondary could 
not be accreted onto the white dwarf and would rather be 
ejected (as noticed by P98). WZ Sge could then be similar 
to another DQ Her star, AE Aqr which is a discless 'ejec- 
tor' (Wynn, King & Home 1997). Such a solution would be 
rather attractive: two, of three known, DQ Her stars would 
be discless magnetic ejectors; apparent similarities between 
WZ Sge and AE Aqr have already been already pointed out 
in P80. 

We think, however, that such a model cannot describe 
WZ Sge. First, Wynn et al. (1997) show that in AE Aqr 
the Hq, Doppler maps of the system are consistent with 
their simulations of a discless flow. Similar maps for WZ Sge 
(Spruit & Rutten 1998) look completely different and clearly 
show the presence of a disc and not just a ring, although 
the authors point out that the brightness distribution is less 
'disclike' than in other CVs. One should bear i n min d, how- 
ever, that as pointed out by Spruit & Rutten (1998) the H a 
brightness may not be a good tracer of matter distribution. 
On the other hand Mennickent & Arenas (1998) find that 
the accretion flow in WZ Sge forms a 'ring' with a ratio of 
the inner to the outer radii ~ 0.3. 




Mass accretion rate (g/s) 

Figure 2. Shown are the critical radius (Eq. 6, with a = 0.01; 
dash-dotted line) and the magnetospheric radius (Eq. 4) as a 
function of the mass accretion rate for a white dwarf mass 
A^wd = 0.45 Mq. The magnetospheric radius relations are nor- 
malized so that the magnetospheric radii in quiescence are equal 
to the Kepler radii which are derived from the observed 28.96 s 
and 29.69 s pulsations. Also indicated are the white dwarf radius 
(ij w( j; dashed line) and outer disc radius (Rdisci dash-dot-dotted 
line) as derived by Smak (1993) and the corotation (Rn) and 
Kepler radii at which the beat between the white dwarf spin fre- 
quency and the Kepler frequency gives rise to the observed os- 
cillations at 28.96 s and 29.69s (dashed lines). The uncertainties 
in i? w< j and -Rdisc are given by the grey hatched areas. The mass 
accretion rates of WZ Sge during quiescence (~210 15 gs _1 ) and 
during the peak of its 1978 outburst (~3 10 18 gs~ x ) as estimated 
by Smak (1993) are indicated with dotted lines. 
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Figure 3. As Figure 2, but now for M wt i = 0.8 Mq. i? w( j is 
determined from the mass-radius relation appropriate for helium 
white dwarfs (Nauenberg 1972). No solution of -Rdisc exists in the 
literature for M wd = 0.8 M . 
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Mass accretion rate (g/s) 

Figure 4. As Figure 2, but now for M w< j = 1.2 Mq. R w< j is 
determined from the mass-radius relation appropriate for helium 
white dwarfs (Nauenberg 1972). -Rdisc (dash-dot-dotted line) is 
the solution given by Spruit Sz Rutten (1998); its uncertainty is 
indicated by the grey hatched area. 



Second, contrary to AE Aqr, WZ Sge is a dwarf nova so 
that we expect an accretion disc to be present prior to, and 
during the outburst (this is independent of the validity of 
the DIM - there is ample observational evidence that dwarf 
nova outbursts require the presence of accretion discs around 
white dwarfs). One could imagine that enhanced mass trans- 
fer prior to outburst could squeeze the magnetosphere and 
allow the reappearance of an accretion disc. However, even 
if this happened, the superoutbursts observed in WZ Sge 
would not occur. This is shown in Figures 2-4 where we 
have plotted the magnetospheric radius as a function of the 
accretion rate for different primary masses. For the magne- 
tospheric radius we use the approximate formula (see e.g. 
Frank et al. 1992) 



R* 



9.8 x 10 8 M~^ /7 M~ 1/7 fi% 7 cm 



(4) 



where M15 is the accretion rate in 10 15 g s _1 and /1 = 
^30 (10 30 gauss cm 3 ) = B,R\ is the magnetic moment of 
a white dwarf with a surface field £>,. Magnetic moments 
corresponding to the 29.69 s oscillation are: ^130 = 34 for 
Mi = 0.45, Mso = 55 for Mi = 0.80 and ^30 = 77 for 
Mi = 1.2 (see also Section 3.4). 

Near maximum of the 1978/1979 outburst most of the 
disc material was accreted onto the white dwarf. Accord- 
ing to Smak (1993) the maximum accretion rate during the 
1978 outburst was 3.2 x 10 18 g s -1 . For primary masses 
larger than Smak's Mi = 0.45, this maximum would be 
lower (Smak - private communication). As the mass accre- 
tion increases during the rise to outburst the magnetospheric 
radius moves inward. Accretion occurs only when this radius 
moves inside the corotation radius. Since no clear oscillations 
are seen during most of the outburst (Patterson et al., 1981) 
this suggests that in WZ Sge the disc was able to reach the 
white dwarf surface (P98). It can be seen, however, from 
Figures 2-4 that at outburst maximum the magnetospheric 



radius given by Eq. (j2J), i.e. the inner disc radius, would just 
reach the corotation radius. For a given maximum accretion 
rate this is the case for all white dwarf masses (because the 
relevant radii are Keplerian) but in reality, for higher masses, 
the magnetospheric radius would not even get there, since 
in this case the maximum accretion rate would be smaller 
than the value plotted on Fig. 2-4. It is obviously impos- 
sible to have the maximum of accretion luminosity at the 
very moment at which accretion onto the white dwarf just 
begins. 

We conclude therefore that WZ Sge is not in a discless 
ejector phase. As we will show below, WZ Sge is most prob- 
ably in a state intermediate between those of the other two 
DQ Her stars: AE Aqr - a pure discless ejector (Wynn et 
al., 1997) and DQ Her itself w hich seems to have a steady 
accretion disc (Warner, 1995a). 



3.2 Blobs at the outer disc edge 

As can been seen in Figure 2, the Keplerian radius corre- 
sponding to the frequency giving the 28.96 s beat with the 
white dwarf spin frequency is very close to the outer disc ra- 
dius as determined from modelling observational data. For 
Mi = 0.45 it coincides, within the error bars, with Rd and 
one can expect that this will be true also for white dwarf 
masses smaller than ~ IMq. We propose, therefore, that 
the 28.96 s beat frequency results from reprocessing of the 
spinning, magnetized, white dwarf beam on plasma blobs 
orbiting at, or close to, the outer disc radius. According to 
Spruit & Rutten (1998) the 'tail' of the observed hot spot's 
H a emission is probably due to material orbiting at the outer 
edge of the disc at the Keplerian velocity. This material had 
crossed the stream-disc interaction region undergoing a se- 
quence of heating and cooling events that are likely to form 
a non-uniform flow. In a stationary accretion disc the outer 
radius is well denned and constant in time so that one can 
expect the motion of plasma blobs there to give coherent 
oscillations. Warner (1995b) suggested that QPOs observed 
during some dwarf nova outbursts could be due to blobs or- 
biting at the outer disc radius. During outbursts, however, 
the outer disc radius moves (in and out; see e.g. Smak 1984a) 
and so coherent oscillations should not be expected in such 
a case. 

Our identification of the 28.96 s period origin is there- 
fore qualitatively consistent with the Spruit & Rutten (1998) 
model of WZ Sge. A quantitative agreement will be diffi- 
cult to achieve because the interpretation of the H a 'tail' as 
due to matter in Keplerian motion leads to their high white 
dwarf mass(1.2 Mq). With such a mass the Keplerian pe- 
riod of the outer disc rim cannot be ~ 733 s. The value of 
the white dwarf mass in WZ Sge is, however, subject to con- 
troversy (see e.g. the discussion in Spruit & Rutten, 1998) 
but for our purpose (see below) the value of 1.2 Mq is too 
high. A lower value of ~ 0.8 Mq would be consistent with 
our hypothesis concerning the origin of the 28.96 s period 
and the estimates of P98. 



3.3 A (marginally) stable accretion disc ? 

As shown by Lasota et al. (1995) the very long outburst 
recurrence time of WZ Sge can be reconciled with the 'stan- 
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dard' values of the viscosity parameter a if the disc is trun- 
cated so that between outbursts it is (marginally) stable. As 
we will see below, at present, the flow of transfered matter in 
WZ Sge does not, probably, form a standard accretion disc. 
We know however that in the past there was an accretion 
disc in WZ Sge (and we expect one to form in the future) 
so it is interesting to see what is the inner radius required 
by the stability criterion. 

The critical value of accretion rate Merit below which 
a disc is in cold, stable equilibrium can be expressed as 
(Hameury et al., 1998; for similar expression see e.g. Ludwig, 
Meyer-Hofmeister & Ritter 1994; Smak, 1984b): 



»> in.-'* 1 0-89 r>2.67 

Merit, 15 = 4.0 a M 1 R 10 , 



(5) 



the accretion disc radius in units of 10 cm. Note that Merit 
depends only weakly on the value of a (in similar formulae 
by other authors the critical M is independent of a [e.g. 
Ludwig et al., 1994] - it should be kept in mind that Eq. (pi) 
is a fit to numerical results). 

It follows from Eq. (|^) that for a given M, an accretion 
disc will be (marginally) stable for r iJ i? C rit where 



D n c »>0.375 0.015 

-Kcrit.io = 0.6 M 15 a 



(6) 



The corresponding Kepler frequency at R cr it is given by 

n nor; nV-0.56 -0.022 /n\ 

^crit = 0.025 M 15 a . (7) 

This equation shows that aj C rit is independent of the 
mass of the white dwarf and only very weakly dependent 
on a. Since the quiescent M in WZSge is ~ 10 15 gs _1 
(Smak, 1993) this leads to a Kepler period at Rait, i.e. 
P(Rcrit) ~ 330-370 s for values of a ~ 0.01-1. Therefore, in 
the framework of the truncated disc model, any beat period 
between the white dwarf spin period and a Kepler period of 
matter in the accretion disc in the quiescent state of WZ Sge 
should have periods shorter than ~30.4s, as observed. 

We notice that the recently re-appeared oscillation at 
29.69 s (P98) is very close to this value. We therefore suggest 
that the 29.69 s period is the beat between the white dwarf 
spin period and the Kepler period of matter near J? C rit- Since 
the inner disc radius, Ri n is larger than i? C rit in quiescence, 
the 29.69 s would mark the inner part of the disc, where 
matter is dominated by the magnetic field. 

We propose that the magnetospheric radius is located 

at 



Rm (29.69) = 8.87 x lCTM?" cm 

corresponding to the Keplerian period Pm 
fastness parameter is now 



Pk 
P s 



16.3, 



(8) 

454.65 s. The 



(9) 



still much larger than 1 so that most of the mass transfered 
from the secondary cannot be accreted by the white dwarf, 
but is ejected from the system (see below). 

Not all matter is ejected - a few percent finds its way to 
the white dwarf and accounts for the pulsed emission. The 
configuration we have in mind is different from disc accretors 
or disc accretors with overflowing accretion stream (see e.g. 
Frank, King & Lasota, 1987; Hellier, 1996). It is also different 
from the pure, discless ejector models of King (1993) and 
Wynn & King (1995). It could be a truncated 'excretion 



disc'. It might be a structure one obtains in the Wynn & 
King (1995) model for low values of the drag coefficient fco. 
According to these authors fco oc B 2 P OT h so that fco for WZ 
Sge should indeed be lower than in AE Aqr. 

From Eqs. (^) , (^) and the white dwarf mass-radius re- 
lation one can deduce that for higher (<; O.7M0) white dwarf 
masses, at outburst maximum the magnetosphere would not 
reach the white dwarf surface. Even taking into account the 
rather crude description of the disc - magnetic field interac- 
tion used above, one can conclude that our model is incom- 
patible with high white dwarf masses. 

3.4 The magnetic field of WZ Sge 

If the inner disc radius is given by Ri n = Rm with Rm set 
by the observed 29.69 s oscillation, we derive magnetic field 
strengths of ~ 3.6 x 10 4 G, ~ 1.6 x 10 5 G, and ~ 1.3 x 10 6 G, 
respectively for Mi = 0.45, Mi = 0.8 and Mi = 1.2. In the 
last case the white dwarf in WZ Sge would have a magnetic 
field surface intensity comparable to that of most IPs. The 
magnetic moment of the white dwarf in WZ Sge is therefore 
3 - 8 x 10 31 G cm 3 . 



fj, ~ 

If the magnetospheric radius is given by Eq. 
for all white dwarf masses of interest, 

Rm > 4J?i 



then, 



(10) 



which, according to Wickramasinghe, Wu & Ferrario (1991) 



(see also Warner, 1995b) means that WZ Sge should be a 
hard X-ray emitter, as it is indeed observed. In this respect, 
as mentioned before, WZ Sge is different from the other two 
DQ Her stars, but for a reason that is obvious: a much lower 
mass transfer rate. 

On the other hand, in outburst, WZ Sge is not a DQ Her 
star but enters the dwarf nova oscillation (DNO) regime in 
which the 'slippage of the surface', i.e. differential rotation 
white dwarf's outer layers is expected (Warner, 1995a,b). 
For WZ Sge the minimum field required to maintain rigid 
rotation is B min ~2x 10 6 M 1 5/2 G, so in WZ Sge B < B min 
(Figure 5). No confirmed DNOs have been see n, however, 
during the last outburst (Patterson et al., 1981). 



4 THE ACCRETION HISTORY OF WZ Sge 

If WZ Sge is in an ejector phase, its white dwarf should 
be spinning down. Indeed, Patterson (1980) reports P s = 
8 x 10" 12 . This value is very large compared with the 
P s — 5.64 x 10~ 14 observed in AE Aqr (de Jager et al., 
1994). The implied spindown power of WZSge would be 
<; 5 x 10 35 erg s _1 , two orders of magnitude higher than 
in AE Aqr. Wynn et al. (1997) find that (assuming that 
plasma interacts with magnetic field in the form of blobs) 
L S pi n scales with M. Since in WZ Sge the mass transfer rate 
is two orders of magnitude lower than in AE Aqr the two 
orders of magnitude higher spindown power is hard to un- 
derstand. One should note, however, an interesting point: in 
WZ Sge the light cylinder radius R c = c/uo — 1.3 x 10 cm is 
larger than the size of the system (a ~ 3.4 x 10 10 cm) which 
allows, in principle, acceleration of particles outside the sys- 
tem. It may well be that results for systems in which R c ^ a 
do not apply to WZ Sge. For particles accelerated to speeds 
close to the speed of light 
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White dwarf mass (solar mass) 

Figure 5. The magnetic field (continuous line) of WZ Sge calcu- 
late assuming that the magnetospheric radius is given by Eq. (0) . 
Nauenberg's mass-radius relation is also assumed. The dashed line 
represent the minimum magnetic field B m i n required to maintain 
rigid rotation of the white-dwarf outer layers during outburst. 



MR 2 c uj 2 « 9 x 10 35 Mis erg s" 1 



(11) 



so that to obtain the spindown power that corresponds to 
P ~ 8 x 10" 12 the escape velocity of particles removing 
rotational energy would have to be a significant (~ 90%) 
fraction of the speed of light. This is, probably, not totally 
absurd since TeV 7-rays have been observed from AE Aqr 
(but not from WZ Sge which is closer) but we would still 
expect a spin-down rate in WZ Sge closer to (or lower) than 
that measured in AE Aqr. A reliable P would be an impor- 
tant test of any model assuming that WZ Sge contains a fast 
magnetic rotator. 

In any case we have to consider how the white dwarf 
in WZ Sge had been spun up to its present rate. The only 
way to achieve this is through accretion disc spin-up. The 
white dwarf spin rate increases as it accretes the Keplerian 
angular momentum of matter at the inner disc edge. An 
equilibrium period is reached when angular momentum is 
accreted at the same rate as it is centrifugally expelled by 
the spinning white dwarf (a magnetic field is not necessary 
for this to work). In the case of a magnetized white dwarf 
the equilibrium period is equal to the Keplerian period of 
the magnetosphere (lu s — 1). The equilibrium spin is then 



P eq = 360Mj~ ' M^j V31 s (12) 

We see immediately that an accretion rate ~ 10 17 g 
s _1 is required to bring WZ Sge to its present rapid rotation 
rate. Such accretion rates are typical of nova-like CVs at 
orbital periods longer than 3 hours whereas WZ Sge is very 
close to the minimum period for CVs where secular mass- 
transfer rates are two orders of magnitude lower. The present 
mass transfer rate in WZ Sge is in good agreement with that 
predicted by mass loss from the secondary being driven by 
gravitational radiation alone. It might even be possible that 
WZ Sge is already on the 'other side' of the minimum period 
(see e.g. Patterson, 1998) . At orbital periods close to 80 min- 
utes the secondary, which due to mass-loss is out of thermal 
equilibrium, stops to contract in response to mass loss and 



the binary system starts expanding to longer orbital periods. 
This 'bounce' is helped by the secondary becoming degener- 
ate (see e.g. King, 1988). It would seem, therefore, that high 
accretion rates are impossible at short orbital periods. 

The situation, however, is not as hopeless as it would 
seem. We are interested here in time-scales of the order of the 
white dwarf spin-up time, i.e. 10 4 — 10 s years. Fluctuations 
on such short time-scale do not modify the secular evolution 
of the binary. 

There is at least one piece of evidence which shows that 
the accretion history of WZ Sge might have had episodes 
of high mass-transfer. There exist, at orbital periods be- 
tween 79 and 92 minutes systems which are supposed to 
have high (~ 10 17 g s _1 ) mass transfer rates. These are ER 
UMa systems which show extremely short intervals between 
superoutbursts (19 - 44 d) and very short (3 - 4 d) 'nor- 
mal' outburst interval. They seem to be the high accretion 
rate equivalent of the low accretion rate systems showing 
superoutbursts: the SU UMa's (WZ Sge is a SU UMa sys- 
tem showing superoutbursts only). The idea that DI UMa 
(one of the 4 known ER UMa's with P or b = 79 min), which 
is more luminous than WZ Sge by a factor ~ 50, spends 
most of its life as 'an ordinary WZ Sge star', but was caught 
during an 'upward surge in accretion', has been proposed 
by Patterson (1998). He speculates that the short period 
CV mass-transfer cycles in which M surges to high values 
(~ 1 — 5% of the cycle time) are somehow associated with 
classical nova eruptions because in addition to the 4 ER 
UMa's, 4 other systems are 'too bright' for their orbital pe- 
riods and 3 of them are classical nova remnants. 

Whatever the reason, what is important for our argu- 
ment is the presence of systems which accrete at high rates 
but have orbital periods close to that of WZ Sge. With such 
accretion rates the white dwarf in WZ Sge can be spun up 
to P cq = 27.87 s. 

When the white dwarf gets to this equilibrium spin rate 
the mass transfer is stable (Wynn & King, 1995). After 
some time the mass transfer rate will return to its secu- 
lar value. This value, however, is two orders of magnitude 
lower than that required for P oq ~ few x 10 s. In a very 
short time (viscous time of the disc ~days) the magneto- 
sphere will start expanding and becoming larger than the 
corotation radius. The system enters into the ejector phase. 
As shown by Wynn & King (1995) this phase is dynami- 
cally stable (for 0.05 Ss q Ss 1). It will stay in this, WZ 
Sge-phase (for about 10 5 years) until it gets to a new spin 
equilibrium corresponding to the low accretion rate (~ few 
minutes). Such a system has been recent ly observed: IP RX 
J0757.0+6306 (Tovmassian et al., 1998) has an orbital pe- 
riod very close to that of WZ Sge (81 ± 5min) but its spin 
period 8.52 ± 0.15 min corresponds to the equilibrium value 
as given by Eq dl2] ) for the expected secular mass-transfer 
rate. The white dwarf's magnetic field in this system could 
have, therefore, a strength close to that of WZ Sge. The 
binary RX J0757.0+6306, therefore, could be a WZ Sge- 
type system at a particular phase of its spin-up/spin-down 
history. (Note that this system is different from the other 
IP with a very short orbital peri od (84 min), RX J 123 8-38 



There the spin period is 36 min (Buckley et al., 1998) so it 
would be a 'classical' IP with fi^i ~ 10.) 

A jump in the mass transfer rate could bring such 
an equilibrium system back to a spin-up phase. If RX 
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J0757.0+6306 has indeed a magnetic field close to the one 
we obtained for WZ Sge and if its mass-transfer rate is close 
the secular one, its accretion disc should be truncated at 
a radius close to the value given by Eq. (^) and therefore 
stable with respect to the thermal-viscous instability. It is 
not yet clear if this system shows any kind of dwarf-nova 
activity. 



5 SUMMARY OF THE MODEL 

We can summarize our model (or 'scenario') as follows: WZ 
Sge contains a 0.8Mq magnetized white dwarf spinning 
with a period of 27.87 s. The magnetic moment is ^ 5 x 10 31 
G cm 3 , or magnetic field at the surface $2 x 10 5 G. 

In quiescence the mass transferred from the secondary 
forms a disc-like structure with an outer radius and inner 
radii at 1.22 and 8.87 x 10 9 M 1 1/3 cm respectively. The accre- 
tion flow is disrupted by the rapidly rotating magnetosphere 
and since the fastness parameter is 3> 1 (the corotation ra- 
dius is much smaller than the magnetospheric radius) most 
of the matter cannot be accreted by the white dwarf. A small 
fraction (Ss 10%) finds its way to the white dwarf surface 
where it is responsible for the observed X-ray emission. 

The white dwarf is in a spin-down phase and most of the 
transferred matter is ejected from the system. The disc-like 
structure formed by the flow might be composed of plasma 
'blobs'. At outburst, which is triggered by an enhanced mass 
transfer event, an accretion disc is (re)created. At outburst 
maximum the disc reaches down to the white dwarf's sur- 
face. 

The white dwarf in WZ Sge will spin down to an equi- 
librium rotation period of a few minutes. In this phase an 
increase of accretion rate will transform it into an ER UMa 
star. After a spin-up phase the system will become a DQ 
Her star once more. 
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